Disordered electron liquid in double quantum well heterostructures: 
Renormalization group analysis and dephasing rate 
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We report a detailed study of the influence of the electron-electron interaction on physical observ- 
ables (conductance, etc.) of a disordered electron liquid in double quantum well heterostructure. 
We find that even in the CEise of common elastic scattering off electrons in both quantum wells, the 
asymmetry in the electron-electron interaction across and within quantum wells decouples them at 
low temperatures. Our results are in quantitative agreement with recent transport experiments on 
the gated double quantum well AlajGai-^^As/GaAs/AIa^Gai-a: As heterostructures. 
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I. INTRODUCTION 

Disordered two-dimensional (2D) electronic systems 
have been remaining in the focus of experimental and 
theoretical research for more than three decades ii The 
experimental discoverjj2i^ of the metal-insulator tran- 
sition (MIT) in a high mobility silicon metal-oxide- 
semiconductor field-effect transistor (Si-MOSFET) in 
1994 became a challenge to a theory. Although during 
last decade the behavior of resistivity similar to that of 
Ref. [HQ has been found experimentally in a wide va- 
riety of two-dimensional electron systems?^ the MIT in 
two dimensions still calls for deeper theoretical and ex- 
perimental understanding. 

Very likely, the most promising theoretical framework 
for studying the 2D MIT is provided by the effective low- 
energy theory, initially developed by Finkelstein, that 
combines the diffusive dynamics due to disorder and 
strong electron-electron interactioui^ Moreover, it is the 
Finkelstein theory that suggested metallic behavior at 
low temperatures long before the experimental discov- 
ery of the MIT in a Si-MOSFET,^i^ Recently, Punnoose 
and Finkelsteinii have shown a possibility for the exis- 
tence of the MIT in the special model of 2D electron 
system with SU{J\f) degrees of freedom in the limit of 
the large number of multiplets. A/" — >■ oo. On the other 
hand, the current theoretical results^i^ do not support 
the existence of MIT for electrons interacting in the sin- 
glet channel only {Af = 1). Therefore, the presence 
of additional degrees of freedom (spins, valley isospins 
etc.) plays a crucial role for the existence of the MIT 
in 2D disordered electron systems. In fact, the impor- 
tance of the multiplet channels of the interaction has been 
confirmed experimentally in Si-MOSFET where a weak 
magnetic field applied parallel to the 2D plane changes 
the behavior of resistivity from metallic to insulating at 
low temperaturesj^^— These experimental findings have 
been explained in framework of the Finkelstein theory in 
the presence of Zeeman and valley splitting The effect 
of intervalley scattering has been taken into account as 



Recently, the Finkelstein theory for disordered elec- 
tron liquid in Si-MOSFET has been subjected to a de- 
tailed experimental check. In particular, the metallic 
behavior of resistivity not far away from the MIT^i^ 
the increase of interaction parameter in the multiplet 
channelsr^"— and the two-parameter scaling near MIT- 
have been observed in experiments. Such the analysis in 
Si-MOSFET is complicated by the presence of (uncon- 
trolled) large valley splitting and intervalley scattering 
rate, A„ « 1/t„ « IK^^ 

As known^^i in n-AlAs quantum wells, 2D electrons 
can also populate two valleys. In addition to Si- 
MOSFETs this system offers opportunity for an exper- 
imental investigation of the interplay between the spin 
and valley degrees of freedom. Using a symmetry break- 
ing strain to tune the valley occupation of the 2D electron 
system in the n-AlAs quantum well, as well as a parallel 
magnetic field to adjust the spin polarization, the spin 
- valley interplay has been experimentally studiedi^i24 
However, the electron concentrations in the experiment 
were at least three times larger than the critical one cor- 
responding to the MITi^ Therefore, the spin-valley in- 
terplay in n-AlAs quantum well has been studied only in 
the region of a good metal, very far from the MIT. 

Disordered electron liquid in double quantum well het- 
erostructures represents a 2D system in which electrons 
in addition to spin have the other degree of freedom: the 
isospin associated with a quantum well. In spite of a 
number of interesting physical phenomena observed in 
electron liquids in double quantum well heterostructures 
without and under strong magnetic field, e.g.. Coulomb 
dragr^ Bose-Einstein condensation of excitona^^, ferro- 
magnetio^i and canted antiferromagnetic phases;^ the 
metal-insulator transition has not been yet addressed ex- 
perimentally. Transport of electrons in double quantum 
well heterostuctures has been studied experimentall y^^ i 
only in the metallic regime far from the region in which 
MIT is expected. 

Recently, detailed experimental research on the in- 
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terference and interaction corrections to conductance 
of electrons in a double quantum well heterostructure 
has been performed i^i^ In particular, two very distinct 
physical situations have been investigated: i) both quan- 
tum wells have equal electron concentrations and mobil- 
ities; ii) one quantum well remains with almost the same 
electron concentration as in case i), whereas the other is 
empty by applying the gate voltage. Surprisingly, it was 
found that the dephasing rate and interaction correction 
to the conductance are almost the same for cases i) and 
ii). 

In the present paper, motivated by the experiments 
of Rcfs. Sllls^ . we develop the theory of the disor- 



dered electron liquid formed in a heterostructure with 
two almost identical quantum wells. We concentrate on 
the case of equal electron concentrations and mobilities 
in both quantum wells [corresponding to the case i) of 
Refs. lslUs^ l. This case will be termed as balance in what 
follows. 

We restrict our study to temperatures (T) satisfying 

the following condition: l/r_| ,As, As as ^ T 1/Tti-. 

Here l/r_| stands for the rate of elastic scattering be- 
tween symmetric and antisymmetric states in the dou- 
ble quantum well structure, Asyis the splitting of these 
symmetric and antisymmetric states, the Zecman 
splitting, and rtr the clastic transport mean free time. 
The temperature behavior of the interaction correction 
to the total conductance is governed by one singlet and 
15 multiplet diffusive modes. We find that the latter 
splits into three inequivalent groups of one, six, and 
eight modes. This grouping occurs due to asymmetry in 
electron-electron interactions across and within quantum 
wells which breaks the rotational symmetry in the com- 
bined spin and isospin spaces [SU (4)]. This reduced sym- 
metry is a distinctive feature of double quantum well het- 
erostructures at the balance and is absent in two-valley 
systems in Si-MOSFETs and n-AlAs quantum wells. We 
identify all relevant interaction parameters and estimate 
their dependence on the distance between the quantum 
wells. To describe the system at low temperatures and 
beyond interaction corrections to conductance, we derive 
the non-linear sigma model and study its renormalization 
in the one-loop approximation. As we demonstrate, the 
renormalization group equations describing the length 
scale dependence of the total conductance and interac- 
tion parameters drive the system towards the fixed point 
corresponding to two separate quatum wells. In spite of 
the symmetry breaking between 15 multiplet modes, the 
renormalization group equations predict the metallic be- 
havior of the conductance at low temperatures. Finally, 
we generalize the expression for the dephasing rate of 
electrons due to the presence of electron-electron inter- 
action known^i^ for a single quantum well to the case 
of double quantum well heterostructures. We find that 
our results are in good quantitative agreement with ex- 
perimental data of Refs. [siUs^ . 

The paper is organized as follows. In Section |TT] we 
introduce the microscopic Hamiltonian, identify relevant 



interaction parameters, study its dependence on the dis- 
tance between quantum wells and introduce the nonlin- 
ear sigma model that describes the low-energy excitations 
in the disordered interacting electron system. Then, in 
Sec, mil we consider the renormalization of the nonlinear 
sigma model in the one-loop approximation, derive cor- 
responding renormalization group equations, and discuss 
renormalization group flow. We derive expressions for 
the dephasing rate due to electron-electron interactions 
in Sec. IIVI Next in Sec.|V]we perform detailed compari- 
son between our theory and recent experimental data on 
transport in double quantum well heterostructures. We 
end the paper with conclusions fSec. IVI|) . 



II. FORMALISM 

A. Microscopic Hamiltonian 

We consider 2D interacting electrons in double quan- 
tum well heterostructures in the presence of quenched 
disorder at low temperatures T <g; r^"^. In the case of 
two almost identical quantum wells an electron annihila- 
tion operator can be written as a linear combination of 
symmetric and antisymmetric states: 

Here electron motion along z axis is confined by the quan- 
tum wells, r denotes a vector in plane perpendicular to 
the z axis, and R ~ r + zBz. The superscript a ~ ± 
denotes electron spin projection, t = ± enumerates sym- 
metric and antisymmetric (— ) states in the double 
quantum well structure and "0? is the annihilation oper- 
ator of an electron with the spin and isospin projections 
equal cr/2 and r/2, respectively. The normalized enve- 
lope function ipi^r{z) — (p(z ± d/2) corresponds to the 
wave function of an electron localized in a single left / right 
well. In what follows, we assume a negligible overlap 
between the states in two quantum wells: the width of 
an electron state in a quantum well [J dz ip'^{z)]~^ ^ d 
where d is the distance between the centers of the quan- 
tum wells. 

In the path-integral formulation, interacting electrons 
in the presence of the random potential V{R) are de- 
scribed by the following grand partition function 



Z : 



(2) 



with the imaginary time action (/? = 1/T) 



S^- j^dt^j dr^'',{rt)[dt+no]rr{rt) - Cai, - Ci^t). 

(3) 

The single-particle Hamiltonian 

V2 1 

Ho = -T /i+-(A,(j + AsAsT) (4) 

znie 2 
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describes a 2D quasiparticle with mass nie- Magnetic 
field B perpendicular to the z axis induces the Zee- 
man splitting As — gfisB. The energy difference be- 
tween symmetric and antisymmetric states in a dou- 
ble quantum well structure yields the splitting Ag^g ~ 
2ip{d/2)ifi' {d/2) /rrie^ The chemical potential is denoted 
as /i, g stands for the effective electron g-factor and /is 
the Bohr magneton. The single-particle Hamiltonian 
is completely analogous to one for a Si(001)-MOSFET. In 
latter case, index r enumerates valleys and Ag^s plays 
a role of a valley splitting. 
Next, the term 



dis 



drr,(rt)Vr,r,{r)rrArt) 



(5) 



describes electron scattering off a random potential 
V{R). It involves matrix elements 



(6) 



In general, the matrix elements Vt-^t2 induce transitions 
between symmetric and antisymmetric states in a double 
quantum well structure. In the case of symmetric random 
potential: V{r,z) = V{r,—z), the system is protected 
from the symmetric-antisymmetric scattering. 

In accordance with the experimental conditions re- 
ported in Ref. (3lll32| , we assume that impurities are con- 
centrated in the middle between two quantum wells. We 
suppose that the random potential created by impurities 
has the Gaussian distribution, and 

{V{R)) ^ 0, {V{Ri)V{R2)} = W{\r, - r^l \z,\, |z2|), 

(7) 

where W decays as the function of its variables at a typ- 
ical distance dw- If the condition 



dz ip'^{z) 



< d, 



(8) 



holds, we can neglect the small difference [proportional 
to ip{d/2)ip' {d/2)], between symmetric-symmetric and 
antisymmetric- antisymmetric scattering rates. Then 

{Vr,r.{ri)Vr,rAr2)) ^ W {\ri - r^l d/2, d/2) Sr.rjr.r,- 

(9) 

Provided correlations in W are short-ranged;^ we find 

1 



(K,..(ri)K3r.(r-2)> = 
1 



5riT2ST3TiS{ri ~ r2), (10) 



2TTl'Ti 

= 27rzy / d^rW{\r\,d/2,d/2). 



Here ly is the thermodynamic density of states of 2D elec- 
trons (including spin). We emphasize that electrons in 
both quantum wells are subjected to correlated disorder 
since they scatter off the very same random potential. 
Recently, under such assumptions, the transconductance 



of a double quantum well structure (the Coulomb drag 
effect with correlated disorder) has been studied by one 
of the authorsi^ 

The small asymmetry in the impurity distribution 
along z axis will lead to the scattering between sym- 
metric and antisymmetric states in the double quantum 

well structure. Its rate can be estimated as 1/t_| ^ 

{b/d)'^/Ti 1/Ti where 6 is a typical length characteriz- 
ing asymmetry. We neglect l/r_| in what follows. 

The interaction part of the action Q reads 



i J dRdR'p[Rt)U{\R-R'\)p[R't) (11) 



where U{R) ~ e^/eR. The dielectric constant is de- 
noted as e. Expanding the density operator p{Rt) = 
{rt)tp'^^ {rt)(pri {z)ipT2 i^) and assuming again that con- 
dition (O holds we obtain 
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/ drdr'^^lirt)ij^^{rt)^P^'^{r't)^P^^^{r't) 

{l + nT2T3n)Un{\r-r'\) 
+ {tiT2 + T3T4)Ui2{\r - r' 



Here 



t/ii(r) = 



dzdz' - 



+ {z- z'Y 



(12) 



(13) 



is the standard Coulomb interaction between electrons 
in a single well. The interaction between electrons in 
different quantum wells 



+ (z - z'Y eVr2 + d"^ 

(14) 

takes into account that electrons are separated by the 
distance d. Due to the difference between Un and U12 
the interaction Lagrangian £int is not invariant under 
global SU{A) rotations of the electron operator t/;^ in the 
combined spin-isospin space. It is the interaction part of 
the action ([3]) that distinguishes the disordered electron 
liquid in double quantum well heterostructures from the 
one in a Si(001)-MOSFET. 

As usual, we single out regions in the momentum space 
of small momentum transfe r^i^i^^i'^^ . Then the low energy 
part of £int can be written as 

Ant = ^ /' ^ ^ab{q)m^\q)m^''{-q), (15) 



Aiy J (27r)2 
m°''(q) = 



a, 6=0 

dk 

(2^ 



iP{k + q)tabrpik). (16) 



Here = {ipX,ip+,^pl,4>_}, tp = 

the 'prime' at the integral sign denotes the integration 
region q {I is the elastic mean free path), and 16 

matrices tab = Ta^cTb stand for the generators of SU{A). 
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Pauli matrices Tq, a = 0, 1, 2, 3 act in the isospin space of 
two wells and Pauli matrices <Tf,, 6 = 0, 1, 2, 3 act in the 
spin space. The matrix of interaction parameters reads 



¥{q) 



fFj Ft Ft Ft\ 

Fs Ft Ft Ft 

F F F F 

yFy Fy Fy Fy J 



(17) 



where 

Ft 



'':(f/n (0))f5, Fy = -Uuil^mFs. 



2 ' ' " 2 
F, ^ u[Uii{q) + Ui2{q)] + Ft, 

Fs^u[Un{q)-Ui2{q)]+Ft. 



(18) 



Here Uii{q) = 27reVge, Ui2{q) = Un{q) cxp{-qd). The 
quantities Ft and Fy are analogous to the standard 
Fermi liquid interaction parameters in the triplet chan- 
nel. They involve averaging of the static part of dynami- 
cally screened interaction U^y^^il^ ^) ^^^^ the Fermi sur- 
face. In the case of equal electron concentrations and 
mobilities in both quantum wells 



(C/n/i2(0)) F5 = 



2lT 



de 

2^' 



UlY/i2i^kFsm{0/2),O) (19) 



where kp is Fermi momentum for a single quantum well. 
The interaction parameter Fs involves the long-range 
part of the Coulomb interaction. In the limit q — >■ it be- 
comes Fs{q) ~ 2k/ q — oo where k = 2TTe^u/e. Within 
the same accuracy, we find 



Fs = >cd + Ft. 



(20) 



At d = — when both quantum wells coincide — the in- 
teraction parameters are equal: Fs = Ft = Fy. Then 
the matrix F corresponds to the case of electron liq- 
uid with two valleys as it occurs in Si(001)-MOSFET. 
In the absence of As^s and the action ([3]) becomes 
invariant under global SU (4) rotations of the fermionic 
fields. In the opposite case of d — >■ oo, the double quan- 
tum well heterostructure is equivalent to two indepen- 
dent single quantum wells. Then we obtain Fs = Fg, and 
Fy = 0. The action (for Asas = A^ = 0) becomes 
invariant under global SU{2) rotations of electron spin 
in each quantum well independently. For intermediate 
values of d, the action ([3]) is also invariant under global 
SU{2) X SU{2) rotations provided Asas and A^ vanish. 



B. Dynamically screened Coulomb interaction 



a b a b a a 



_i_ 3 a 




FIG. 1: Dyson equation for the screened electron-electron in- 
teraction in RPA. Thick wavy line denotes screened interac- 
tion, thin wavy line is bare interaction, solid lines are electron 
Green's functions, and dashed lines are impurity lines. Indices 
a and b can be 1 or 2. Index a equals 1 (2) if index o is 2 (1). 



ically screened interactions 

Uu+Il2[U?i-U?2] 



7-rscr _ 
•^11 — 



TP 



1 + [Hi -f U2]Uii + Uin2[Uf^ - ul^] ' 

1 + [Hi + n2]c/ii + nin2[t/2i - ul^] ' 
1 + [Hi + n2]c/ii + iiiii2[ul^ - u?2] ' 



(21) 
(22) 
(23) 



The polarization operators can be written in diffusive 
approximation as 



Uj{q,uj) 



Djq^ — iu) ' 



J = 1,2 



(24) 



where Dj is the diffusion coefficient in the j-th quantum 
well. We mention that for Di ^ D2 the dynamically 
screened Coulomb interaction in the first well Ul'^{q,uj) 
does not coincide with the one (C/|"((7, a;)) in the second 
well. 

If the electron concentrations and mobilities in the 
quantum wells are the same then Di = Z?2- In this case 

rrsc 
•^11 



rscr 
11 



Vl^ and 



rrscr 
^11 



/12 



2vq\ ^ 
±- 




1 



Dq 



q + x{l- 



(25) 



As one can see, at ^ 1 the effect of the right well 
on the dynamically screened interaction in the left well 
is negligible. In the opposite case, qd ^1 the right well 
affects the dynamically screened interaction in the left 
well only at xd < 1. 



The interaction parameters Ft and Fy involve the 
screened Coulomb interaction. Solving the Dyson equa- 
tions in the random phase approximation (RPA) (see 
Fig. [1]), we obtain the following results for the dynam- 



C. Estimates for the interaction parameters 

Let us estimate the interaction parameters Ft and Fy in 
the case of equal electron concentrations in both quantum 
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wells. By using Eqs. (PT|) and ([^ we find 



Ft ±F,, = - 



^(^l -|- g-2fcFdsine/2'j 



47r 2fcF sin f + ± e-2feF<isi„e/2) ' 



... (2^) 

To justify the RPA which has been used in derivation of 
Eqs. (pijl - ip^ we assume that the condition x/kp <C 1 
holds. As follows from Eq. both Ft and Fy are neg- 
ative and \Ft\ ^ \Fv\- The interaction parameter Fg is 
negative at small d and positive at large d. The depen- 
dence of the critical distance dc at which F^ vanishes on 
the parameter x/kp is shown in Fig. [51 We mention that 
^ \Ft\ ford<4. 
It is instructive to compare the results for Ft, F^ and 
Fg with the case of a single quantum well for which the 
interaction parameter in the triplet channel is given aa^ 



27r 



:ln 



1 + Vl - 



1 - Vl - ^2 



(27) 



In the limit x — )■ the function ^0(2;) acquires the fol- 
lowing asymptotic form 



Go{x) w a;ln(2/x), 



X < 1. 



(28) 



Provided kpd 3> 1, the interaction parameters for the 
case of double quantum wells with equal electron con- 
centrations can be estimated as 



Ft = F^ 



1 



Snkpd 



Gii^d), F,, 



1 



Snkpd 



G2{>id), (29) 



Giix) = + ^ El [ - 



x + 1 

Q2{x) = Gi{x) - 



- 1 

^xe'Ex^x) 



x-\ 



x^\ 



Here E\{x) = dt exp{—t)/t is the exponential inte- 
gral. 

Finally, we mention that the interaction parameters 
Ft and F„ can be estimated (from above) as ^ 
[go{>i/kF)+go{x/2kp)]/{4n) and \F,\ < Goi^^ /kp)/{4n). 
Even for values oi x/kp ^ 1, it yields \Ft\ < 0.3 and 
IF, I 5,0.2. 



D. Non-linear a model 

At low temperatures, Trti- ^ 1, the effective quantum 
theory of 2D disordered interacting electrons described 
by the microscopic action (O is given in terms of the 
non-linear sigma model. The latter describes interaction 
between low-enegy modes which are the so-called "Dif- 
fusons" and "Cooperons" . As well- known , ^'^'^^i'^^ the in- 
terference ("Cooperon") contribution to the conductance 




x/kp 



FIG. 2: Value of the parameter xd at which Fs = versus 



is not sensitive to the presence of As and Ag^s (in the 
absence of ). Furthermore, the interference cor- 
rection is cut off by weak magnetic fields and does not 
influence the scaling of observables with temperature at 
B > l/eDr^. Therefore, we shall ignore the interference 
correction in the intermediate calculations for a sake of 
simplicity and shall discuss its role in Sec. |Vl 

In general, Cooperons are also involved in the interac- 
tion correction to the conductance and the renormaliza- 
tion of other interaction couplings. The corresponding 
contributions are proportional to the interaction param- 
eter in the Cooper channel. For Coulomb interaction, 
the latter is repulsive and remains small in the course of 
the renormalization for 2D electron systems<^ Moreover, 
physically, a moderately weak magnetic field B > T/eD 
applied parallel to the z axis is enough to suppress the 
interaction effects in the Cooper channel 4^ 

Neglecting the Cooper channel, the effective theory in- 
volves unitary matrix field variables Qmr?Tir2'^^ ('') which 
obey the nonlinear constraint Q'^{r) = 1. The integers 
aj ~ 1,2, Nr denote the replica indices. The inte- 
gers m, n correspond to the discrete set of Matsubara 
frequencies e„ = 7rT(2n -I- 1). 

The effective sigma-model action is 

S = S,+Sp+SsB. (30) 
Here Sa represents the free electron parl4^ 



"32" 



dr tr(Vg)^ 



(31) 



with a XX = ^T^v*D denoting the mean- field conductance 
in units of /h. The thermodynamic density of states 
v^, = ra^jix involves an effective mass m» renormalized 
due to interactions. The symbol tr stands for the trace 
over replica, the Matsubara frequencies, spin and isospin 
indices. The Finkelstein term^i^ 

= \ d^Y. ra6trO,,g(r)tr/f„UQ(r-) 

an;ab 

-f-47rrz y dr tr ^{Q - A) - 27rTz ^ dr tr rjA (32) 
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involves the electron-electron interaction amplitudes Tab- 
The bare value of the factor z is determined by the ther- 
modynamic density of states: z — -Kv^j ^. The quantity z 
has been originally introduced by Finkelstcin in order to 
ensure the consistence of the renormalization group equa- 
tions with the particle number conservation)^ Physically, 
the renormalization of z is responsible for renormaliza- 
tion of the specific heat4i and determines the relation 
between the frequency and length scales, thus playing a 
crucial role at the criticality near the MITi^ 

The interaction amplitudes Vah are related to the in- 
teraction parameters Faj, introduced above as^i^i^ Tab = 
— zFab/(l -I- Fab). Therefore, the matrix F has the struc- 
ture similar to the matrix F (see Eq. (flT)) ) and 



Ft 



F„ 



1 + F, 



It 



l + Ft 



l + F, 



(33) 



The matrices A, 77 and are given as 



A"^, = sign {LOn)5nm5°'^tQQ, rf^^^ = nSnmS°''^too, 

{IlrJn = <5„-™,fc<5"^5'^^%. (34) 

The action S„ -\- Sp is invariant under the global ro- 
tations Qj:^irr\r^(r) ^ ul\l-^QZ?A'^,f,\{r)[u-X\l\ with 
u ~ X]a=o X]&=o "^abtab- This rotation correspond to the 
global SU{2) X SU{2) symmetry of the action + Sp- 

The presence of Ag and/or As as generates the sym- 
metry breaking terms. In general, they can be written 



X = n Xn-^n P^^Y ^^e important role4^i^ For exam- 
ple, J^-aigebra allows us to establish the following rela- 
tions 

ii I^tabe'^Qe^'^ = tr/^Ue*^"Qe-'>^« +8m(xafc)"„, 
- 4 ^ n^XabUXabT-,, (36) 

an-^ab 

where xo = X^qXo^o ■ With the help of Eqs. one 
can check that the relation Tg = —z guarantees the so- 
called J^-invariance4^ It is the invariance of the action 
So- + Sp under the global rotation of the matrix Q with 

Xab = X^aoSbO- 



E. Physical observables 

The most significant physical quantities in the the- 
ory containing information on its low-energy dynamics 
are physical observables cr^^, z', and z^^ associated with 
the mean- field parameters axxi z, and Zab of the ac- 
tion (j30p . The observable a'^^ is the total DC conduc- 
tance as obtained from the linear response to an elec- 
tromagnetic field. The observable z' is related with the 
specific heat4i The observables Zab determine the static 
generalized susceptibilities of the 2D electron system^i^ 
as Xab — ^z'abl''^- ^hc couductancc (7^3. can obtained 
from 



SsB = iZab Aab J drtrtabQ + J drAl,. (35) 

For the symmetry breaking by the Zeeman splitting one 
can choose tab = ^03 and A03 = A^. In the case of the 
splitting AsASi the generator tab equals t^Q. Splitting 
Aab set the cut-off for a pole in the diffusion modes ( "dif- 
fusons"). In what follows, we shall be interested in high 
temperatures (T ^ Aab) or, correspondingly, in short 
length scales L ^ ^Dj Aab such that the cut-off is irrel- 
evant and the electron system behaves as if no symmetry 
breaking terms are exist. We shall use the symmetry 
breaking term Ssb only as a source, assuming infinitesi- 
mal Aab. 



1. T -algebra 

The action ([50]) involves the matrices which are for- 
mally defined in the infinite Matsubara frequency space. 
In order to operate with them we have to introduce a cut- 
off for the Matsubara frequencies. One should send the 
cut-off to infinity at the end of all calculations. Then, the 
set of rules which is called J^-algebra can be established.— 
The global rotations of Q with the matrix exp(ix) where 



dr'((tr/;^Q(r)VQ(r)tr/!!„Q(r')Vg(r'))) 



64Dn 



(37) 



after the analytic continuation to the real frequencies: 
iLOn — > w + iO^ at w — > 0. The expectation values are 
defined with respect to the theory ((30|) and D = 2 stands 
for the spatial dimension. The physical observable z' can 
be extracted from the derivative of the thermodynamic 
potential per the unit volume with respect to temper- 
aturei^ 



1 



d n 



27r trryA dTT' 
z'^i^ are given as 



(38) 



2NrdAl, 



(39) 



It is worth mentioning that, alternatively, the observable 
parameters ct^^,, z^^ and z' can be found from the back- 
ground field procedure. 



7 



III. ONE-LOOP RENORMALIZATION 

A. Perturbative expansions 

To define the theory for the perturbative expansions 
we use the "square-root" parameterization: 



Q = W + A\/l - W^, 



W 



w 
wt 



(40) 



The action pop can be written as the infinite series in 
the independent fields w and . At short length scales 
L <C ■\/ <7xx / (zab^ab) which we are interested in, the sym- 
metry breaking term Ssb can be omitted. Then the 
propagators for fields w and can be written in the 
following form 



xDq{uJi2) 



327rTrah 



X<5af,;cd<S"^"^<5"^"^<5„,,,„3,, (41) 

where a;i2 = Em — = 27rTni2 = 2i:T{ni — 71,2) and 



16(z + Tab)(^n 



(42) 



We use the convention that the Matsubara frequency 
indices with odd subscripts 71,1,71,3,... run over non- 
negative integers whereas those with even subscripts 
71,2, ri4, . . . run over negative integers. 



B. Relation of Zab with z and Tab 

The dynamical susceptibility Xab{^, q) which describes 
the linear response of the system to time-dependent sym- 
metry breaking amplitude Aab can be obtained froni^ 

Xab{^^n,q) = ^ - Tzl,{tv I^^tabQ{q)tT rjabQi-q)) 

TT 

(43) 

by the analytic continuation to the real frequencies: 
iujn uj + iO^ . In the tree level approximation Eq. (|43p 
yields 

Xab(«W„, q) = 1 D'^'iujn) . (44) 

The action S^^Sp is invariant under the global rotations 
Q -> uQu^^ with u = I]a=oI^Lo"a6iab- This implies 
that the quantities corresponding to operators ttt,'^'' and 
m}^ conserve, i.e., xo6('^,g = 0) = Xi6(w,q = 0) = 0. 
In order to be consistent with this physical requirement, 
the relations 

Zab = z + r* = z(l +7t), = 0,1,6=1,2,3, 

zw = z + ts= z{l + ^s)- (45) 



should hold. Therefore, renormalization of the interac- 
tion amplitudes T s and Tt can be easily found from, e.g., 
renormalized quantities Zg^ and zJq. However, it is not 
the case for the interaction amplitude = z^y. There 
is no simple relation between r„ and 



Zy = Z2b = Z3b, 



6 = 0, ...,3. 



(46) 



Therefore, the physical observables cr^^, z', 7^, 7^, 
and z'y completely determines the renormalization of the 
theory ([5(1)) at short length scales L <C \/<Jxx/ Zab^ab- 



C. One-loop results 

Evaluation of the conductance according to Eq. ([57)) 
in the one-loop approximation yields 



(^xxii^n) = (Tx 



1287rT 



ab w„,>0 



(27r)i: 



X mm{LL>jn,(^n}Dp{LL>m +UJn)Dp{uJ„i)D^^^\uJm). (47) 

Performing the analytic continuation to the real frequen- 
cies, iuJn — > o; + iO^, one obtains the DC conductance in 
the one-loop approximation: 

ab 

X ^ {n coth ^) [Z)«(f))]^i^f )^«(17). (48) 
Here D^{n) and Dl''''^-^{n) are retarded propagators 



corresponding to Dp{uJn) and D^p''\ijJn), respectively 



- (16(z + Tab)lcrxx) ■ 



(49) 



We mention that the result ([^5)) can be also obtained with 
the help of the background field procedure^ applied to 
the action ([5T]) - ([5^ . 

In order to compute z' , we have to evaluate the ther- 
modynamic potential f2. In the one-loop approximation 
we find 



dT 



z + 



■E 



w„>0 

iz + rab)Dl^''HLOn)~zDp{LOn) 



(27r)i 



Following definition ()38p . we obtain from Eq. ()50p 



ab 



(27r)0 



Dp{Q). 



(50) 



(51) 



Next, we evaluate in the one-loop approximation the gen- 
eralized susceptibility Xab(*w„,q) at g = and a;„ — >■ 0. 



8 



Then, according to Eq. (P^ . we find 



cd;e/ ^ ' 



(52) 



where C^j^.^j denotes the structural constants of SU{A): 
[tcd,tef] = EabCcleftab- Applying Eq. ^ ^OT (ab) = 
(10) and (ab) = (01), and by virtue of relations (|^5|) we 
obtain 



z' + t' ^ Z + ts - 



2^-V(z + r,,)^ 



(27r)i 



(53) 



297r(z + rt)2 /■ d> V- 

x{[l?f )(..„)]% [42°)(..™)]^-2I?^(a.„0}. (54) 

In order to find renornialization of r„, one cannot use 
the static generalized susceptibility since there exists no 
simple relation between Zy and r^,. We use the the 
background-field renornialization procedure (see details 
in Appendix E| and find 



^'ab — Tab 



1 



(27r)» 



cd\ef 



Ld,„>0"' ' cd-ef 



(55) 

Here symbol sp denotes trace over spin and isospin in- 
dices. Using Eq. ([55]) for {ab) = (02) , we obtain 



2i07rr?, 



(27r)i 



(Jxx ^ J (27r)» f 



(wm). (56) 



It is worthwhile to mention that the results ([51]) . ([55)1 
and ([5^ can be also derived from Eq. ([55)1 . Equa- 
tions (gHl) dlSl), (IMl) and ^ allow us to extract 
one- loop rcnormalization of conductance a^x , parameter 
z and interaction amplitudes Eg, Ft and Ft,. 



D. Renormalization group equations 

Applying the minimal subtraction scheme (see, e.g., 
Ref. M) to Eqs. gH), dSil), dMl), (ESI) and we de- 

rive the following one-loop results for the rcnormalization 



group (RG) equations which determine the T = behav- 
ior of the physical observables with changing the length 
scale L in D = 2 dimensions: 



dcjxx 2 
d^s _ 1 + 7. 



[1 + /(%) + 6/(70 +8/(7.)], 

Is - Iv 



TTCTxa; 
d-yt 1 + 7t 

dS, T^CTxx 

djv _ _L 
d^ ~ 



1 - 674 - 7s + 87„ + 167^ 
1 - % + 27t + 87„ 



1+71- 
7t - 7t. " 



din z 



d( 



1 + 7s + 7o - Iv (67* + 7s) + 87^ 
' 7s + 674 + 87„ - 1 



(57) 
(58) 
(59) 
(60) 
(61) 



Here f{x) = 1 - (1 -|- a::~^)ln(l + x), ( ^ \nL/l and we 
omit primes for a brevity. Equations ([ST I - lpO)) constitute 
one of the main results of the present paper and describe 
the system at the length scales L ^ ■\/ Uxx/ (zab^ab)- 

It is worthwhile to mention that the right hand side of 
Eqs. and (|59D is not polynomial in the interaction 
amplitude jy. To the best of our knowledge, the one- loop 
RG equations for interaction amplitudes are quadratic 
polynomials in all cases studied previously] ^' ^'^1 ^^'^^1'^^ 
This fact is deeply related with invariance of the action 
Sa + Sp under the global rotation of the matrix Q with 
the matrix cxp(ix) (see Sec. HID II) . As it follows from 
Eqs. p6p . iSo- + Sp is invariant under such global rota- 
tion with Xab ~ X^acSbd whcrc c = 0, 1 and d ~ 1, 2 or 
3 provided 74 — —1. The same holds for the global ro- 
tation with Xab = X'^ai'^bo if 7s = —1- This invariance 
guarantees that 74 = —1 and 7s = —1 are fixed points of 
the RG equations. Therefore, the latter have to be well- 
defined at 7t = — 1 and 7s = —1. However, for 7„ = — 1 
the action Scr + Sp is not invariant under the global ro- 
tation of the matrix Q with Xab = X^ac^bd with c = 1, 2 
and d = 0,1,2 or 3. It is exactly this noninvariance 
that allows appearance of factors 1/(1 + 71,) (diverging 
at 7„ = -1) in Eqs. ^ and ([5^)). 

The renormalization group equations ([57|) -(|Bn i) pos- 
sess a rich four-dimensional (cra;^;, 7s, 7t, 7i;) flow dia- 
gram. First of all, there is the two-dimensional sur- 
face 7t = 7„ = 7s which is conserved under RG flow. 
It corresponds to the case of coinciding quantum wells 
{d ~ 0). In this case, the RG equations (|57)) - ((6T|) are 
completely equivalent to ones for the two- valley electron 
liquid. However, this two-dimensional surface is unsta- 
ble: a small initial mismatch (e.g., due to finite d) in the 
condition 7t = 7t; = 7s increases during RG flow. Sec- 
ondly, the RG flow conserves the two-dimensional surface 
7„ = 0, 7s = — 1 which is stable. It describes the limit of 
two separate quantum wells (d = 00). In addition, there 
are some interesting features of RG flow. For example, 
there is a two-dimensional surface 7t = 7s = — 1 which is 
conserved by RG flow. There is an accidental fixed line 
7s = —1, 7i, = —1/2, 7t = —1/3. However, these features 
are not accessible in the double quantum well structure. 

Indeed, the initial values of the parameters 7^, 7^ and 
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FIG. 3: Dependence of the parameters 74, 7„, 7s on ^. Initial 
values are 7i = 0.35, 7,, = 0.01, 7s — —0.77, and a^x = 6. 



7s satisfy 



7t > 7u ^0, 7t > 7s. 



(62) 



Then, using Eq. ([55| - (pn|) one can prove that under RG 
flow i) the conditions 7* ^ 7i; ^ and 74 ^ 7s hold; ii) 
"ft always increases. Starting from initial values of the 
parameters 74 , 7„ and 7s satisfying Eq. ((62)) the RG flow 
develops in such a way that 7^, vanishes, 7s tends to —1 
and 7t increases towards infinity as shown in Fig. [31 

The conductance a^x demonstrates metallic behavior 
as in the case of two- valley electron liquid. It increases at 
large length scales. Depending on the sign of the parame- 
ter Kee = l+/(7s)+6/(7t)-l-8/(7^), the conductance can 
develop both monotonic {Kee < 0) and non-monotonic 
behavior {Kee > 0) (see Fig. |4]). The phase diagram for 
the parameter Kee is shown in Fig. [5] At x/kp < 0.4 
the parameter Kee is positive for all values of }<d. With 
increasing x/kp & domain of negative values of Kee de- 
velops at small values of xd. 

The conductance cr^^ defined in Eq. p7p and rcnormal- 
ized in accordance with Eq. (|57p is the total conductance 
of double quantum well structure. In general, one can 
write CT^^ = 0-^1 + (T22 + ^'12 + <^2i, where cr'n and (732 
are the intrawell conductances of left and right quantum 
wells respectively, and a'12 a-nd denote the transcon- 
ductances responsible for a drag effect. At the balance, 
symmetry yields that (t'h = cr'22 and o'5^n = o[^. 

Although in experiments of Refs. (SlIlS^ only the 
total conductivity o'^^ has been measured, such dou- 
ble quantum well heterostructures with correlated dis- 
order at the balance allow for experimental study of 
transconductancc contrary to the two- valley electron sys- 
tem in Si-MOSFET. It was shown2£ that in the presence 
of electron-electron interaction one-loop contribution in 
the particle- hole channel (only "diffusons") to the DC 
transconductancc a'^^ vanishes. As a result, the one-loop 
contribution to the DC transconductancc is entirely de- 
termined by the particle-particle channel ( "Cooperons" ) . 
However, in Ref. [33 only the interwell interactions {U^) 
were taken into account. As shown in Appendix [Bl an 
accurate treatment of both interwell (C/i" ) and intrawell 



(t/^") interactions (i.e. taking into account all interac- 
tion couplings Fs, Fs, Fj, and F„) does not change the 
conclusion of Ref. [s^ : the particle-hole ( "diffuson" ) con- 
tribution to the DC transconductancc vanishes in the 
one-loop approximation. 



IV. DEPHASING RATE 

The presence of the right well changes the properties of 
electrons in the left well. One of the important quantities 
characterizing interacting electrons in a random potential 
is the dephasing rate. Its dependence on temperature de- 
termines the behavior of the weak-localization correction 
to the conductance. In this section, we investigate how 
the presence of the right well changes the dephasing rate 
of electrons in the left well compared to the case when 
the right well is empty. 



Contribution from the interaction in the singlet 
channel 



We start from the case of the interaction in the singlet 
channel only. According to Eq. ([2T|) . electrons in the 
right well screen interaction between electrons in the left 
well and vice versa. The dephasing rate of electrons in 
the left well due to the presence of electrons in the right 
well can be found from the following expression which 
generalizes standard ontS liSl 



1 



dio 



lmUff{q,uj) Diq^ 



TT J (27r)2 sinh(a;/T) 



(63) 

Expression for the dephasing rate of electrons in the right 
well can be obtained from Eq. by substitution of 
U22 and D2 for Uu and Di, respectively. At the balance 
which we are interested in, the dephasing rates in the 
left and right wells arc the same. Under the following 
assumption d, ^ Lp = ^ D jT ^ wc find 



-=-4s- 



In Tr„ 



where is the function of the parameter >cd: 

{udf 



A.. 



1 



(1 + >rd)(2 xd) 



(64) 



(65) 



We mention that in the absence of electrons in the right 
well (formally this case corresponds to the limit d — > 00) 
the dephasing rate is maximal: As = 1. Equation (1651) 
was used for analysis of the experimental data in Ref. [3l| . 



B. Contribution from the interaction in the 
multiplet channels 

In general case, one has to take into account contribu- 
tions to the dephasing rate from the interaction in mul- 
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In the absence of interaction in the multiplet channels, 
i.e., for Ft = Fy = and % = —>cd/{l + nd), this result 
transforms into Eq. (j64p . We mention that the interac- 
tion parameters % , 7* and 7.0 as well as conductance axx 
should be taken at the length scale Lt = y/ (Jxx/ zT. 

It is worthwhile to compare Eq. with the result for 
the dephasing rate in the absence of electrons in the right 
well;^ Taking the limit d — ?> 00, i.e., setting 7^ = — 1, 
7i, = 0, and 74 = 74,0 1 we obtain 



Ho 
2 + 7t,o 



(70) 



FIG. 4: Dependence of the inverse conductance l/a^x on 
^. Solid (red) curve is plotted with initial values 74 — 0.35, 
jv = 0.01, 7s — —0.77, and cjxx = 6. Dashed (blue) curve cor- 
responds to 7t = 0.198, jv = 0.135, 7s — 0.57, and axx = 6. 



where initial value of 74^0 is 7t,o = —Ft/{^ + Pt)- 



COMPARISON WITH THE EXPERIMENT 



0.5 



FIG. 5: The phase diagram for the parameter K^e = 1-1- 
/(7s) -I- 6/(7*) + 8/(71,). It vanishes on the solid (blue) line. 
Kee is negative in the filled region below the solid (blue) curve 
and is positive above. Dashed (red) line indicates Fs — 0. 



tiplet channels 1^ We restrict ourselves to the case of the 
balance. Generalizing the well-known result^ for the sin- 
gle well we can write the dephasing rate in the left well 
as 



1 



where 



(27r)2 sinh(w/r) 



(66) 



U''''^\q,uj) = -^Dl;^^''^iLj)[D^{uj)]~\ (67) 

Performing integration over momentum and frequency, 
we find 



1 T 
— In Tr. 

ZCxx 



with 



It 



7^ 

IV 



-It 



(68) 



(69) 



Recently, the interferenceai and interactional correc- 
tions to the conductivity of the gated double quan- 
tum well Ala;Gai_2,As/GaAs/Ali,Gai_a;As heterostruc- 
tures have been studied. Two heterostructures, 3243 
and 3154, distinguishing by the doping level have been 
investigated. From analysis of positive magnctoconduc- 
tivity the dephasing rate has been extracted. By tuning 
the gate voltage, the electron concentration in the right 
quantum well were controlled in the experiment. 

We consider two characteristic cases: I) electron con- 
centrations and mobilities {M.) of both quantum wells 
are equal: ni = n2 = n and Mi = M2 = -M; II) the left 
quantum well has electron concentration ni = n and mo- 
bility Ml = M whereas the right quantum well has elec- 
tron concentration 712 = 0. The electron concentration n 
has been high such that the conductance (Jxx was about 
80. Therefore, physics described by RG equations ((ST]) - 
(|6ip was not observed. The main unexpected findings of 
Refs. [silfs^ were as follows. Dephasing rates (coefficient 
,4) and interaction correction (parameter K^e) extracted 
in cases I) and II) were practically the same. At first 
glance, it is counterintuitive since there arc 15 multiplets 
in the case I) and only 3 in the case II). 

After Refs. [sills^ l we summarize the experimental val- 
ues of relevant parameters in Table HI The theoretical 
estimates for the interaction parameters in cases I) and 
II) are presented in Table |TT1 As one can see, in the 
experimentally studied case of xd = 3.6 the interaction 
parameter Fy is negligible. Ft and F° coincide with each 
other and Fs is equal approximately to xd. The compar- 
ison between theoretical estimates for K^e, A, ifee,o ^^nd 
^0 with experimental data (whercever it is possible) is 
summarized in Table IIIIl Our theoretical estimates are 
in good quantitative agreement with the experimental 
ones. Our results explain why the interaction corrections 
and dephasing rates in cases I) and II) were found to be 
practically the same in the experimentsi^'^ Since the 
parameter Kee,o is positive for >£'/fcF ^ 1 a drastic effect 
in the interaction correction could be seen by tuning the 
gate voltage from case I) to case II) in double quantum 
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TABLE I: Parameters for samples studied in Refs. [3lll3! 



sample 


#3154 


#3243 


n, 10" cm~^ 


4.5 


7.5 




1.7 


2.2 


X, 10*^ cm^i 


2 


2 


d, IQ-^ cm 


1.8 


1.8 


x:d 


3.6 


3.6 


kpd 


3.06 


3.95 


n/kp 


1.18 


0.91 


Theoretical estimates of interaction 


sample 


#3154 


#3243 


Fs 


3.34 


3.37 


Ft 


-0.26 


-0.23 


F, 


-0.009 


-0.007 


fa 


-0.77 


-0.77 


7t 


0.35 


0.30 


> 


0.009 


0.007 




-0.26 


-0.23 


7t,o 


0.35 


0.30 



well structures with Kd < 1 for which one can expect 
Kee < (see Fig.©. 

As mentioned in the Introduction, our theory is valid 
at temperatures T ^ Asas, '^sA/ta In the experi- 
ments of Refs. [3lll32| the Zeeman splitting (at relevant 
magnetic field which was used in order to extract interac- 
tion correction) and Asas were estimated as As < 0.2K 
and Asyis < IK. A small asymmetry in the impurity 
distribution along z axis presented in the double quan- 
tum well heterostructures used in Refs. SlHs^ leads to 
appearance of scattering rate between symmetric and an- 
tisymmetric states. The corresponding scattering rate 
) can be estimated from temperature and mag- 
netic field dependence of weak-localization (interference) 
correction to conductivity. 



TABLE III: Comparison of theoretical estimates and exper- 
imental findings [Kee,0 = 1 -f 3/(7t,o), A'ee = 1 -f /(7s) 4- 

6/(70 + 8.f (7.)]- 





Theory 


Experiment 




#3154 


#3243 


#3154 


#3243 




0.59 


0.72 


0.50 ± 0.05 


0.57 ±0.05 


Kee.O 


0.52 


0.59 


0.53 ± 0.05 


0.60 ± 0.05 


A 


0.89 


0.86 






Ao 


1.15 


1.12 






A/Ao 


0.77 


0.77 


1.00 ±0.05 


1.00 ±0.05 



As known^i^i^ in the absence of scattering be- 
tween symmetric and antisymmetric states neither A., 
nor A5^5 does not influence the weak-localization con- 
tribution. In the absence of magnetic field, the weak lo- 
calization correction to the conductance in both asymp- 
totic cases Asas *C !/'''+- ^i^d Asas ^ ^ /'''+- can be 
written as 



Sa. 



WL 



1 



In 



if! 



1 

Tl2 



(71) 



where 1/ti2 ~ min{A|^sr^ ,1/t^ }. The tempera- 
ture dependence of the weak-localization correction ([7T|) 
smoothly interpolates between the result known for a 
two- valley system at high temperatures (1/t^ ^ ^I'^vi) 
and the single- valley result at low temperatures <C 
l/''"i2)- In experiments |3l| the characteristic time ri2 
was estimated from the suppression of weak-localization 
correction due to perpendicular magnetic field as l/ri2 < 
Q.IK. Together with the estimate Asas ^ 1^ it implies 
that I/t-i ~ 1/ti2 < O.IK. Therefore our theory is ap- 
plicable at temperatures T > IK. It is this temperature 
range that was studied experimentally in Refs. |3lll32| . 



VI. CONCLUSIONS 

To summarize, we have developed the theory of the 
disordered electron liquid in a double well quantum het- 
erostructure with equal electron concentrations. We have 
identified all relevant interaction parameters and found 
their dependence on the distance between quantum wells. 
To describe the system at low temperatures, we have de- 
rived the interacting non-linear sigma model and studied 
it renormalization in the one-loop approximation. We 
have obtained the renormalization group equations de- 
scribing the length scale dependence of the conductance 
and interaction parameters. We have found that upon 
the renormalization the system fiows towards the fixed 
point corresponding to two separate quatum wells. The 
RG equations predict the metallic behavior of the con- 
ductance. We have evaluated the dephasing rate of elec- 
trons due to the presence of electron-electron interaction. 
This expression takes into account screening of electron- 
electron interaction within one quantum well by electrons 
from the other quantum well. 

We did not consider contributions to the one-loop RG 
equations from the particle-particle ("Cooper") chan- 
nel. The interaction effects related to the Cooper chan- 
nel are governed by the corresponding interaction am- 
plitude which is always small for 2D electron systems 
with Coulomb repulsion, so that one can neglect iti^ 
As for the interference contribution to conductance, for 
1/T(p ^ l/''"i2, it can be taken into account by the sub- 
stitution of 1 + 2 for 1 in the square brackets of Eq. ([57]) . 
This does not change qualitative behavior of the interac- 
tion amplitudes 7s, 7t and discussed above. However, 
the interference contribution makes behavior of the con- 
ductance always non-monotonous. 
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We performed detailed comparison between our theory 
and experimental data^ii^ We explained main experi- 
mental results and found good quantitative agreement. It 
would be an experimental challenge to construct the dou- 
ble quantum well heterostructure with >cd < 1. Then, ac- 
cording to our predictions, one can expect a change from 
non- monotonous to monotonous behavior in conductance 
in the presence of small perpendicular magnetic field (to 
suppress interference contribution) when the right well is 
depopulated by tuning the gate voltage. It would be also 
interesting to experimentally study the Coulomb drag ef- 
fect in such heterostructures with correlated disorder. 

Finally it would be worthwhile to extend our analysis 
to temperatures less than the symmetry breaking energy 

scales Agyis, and l/r_| At such low temperatures 

one may expect different behavior of transport in dou- 
ble quantum well structures as compared to two-valley 
electron systems studied recentlyJ^ii^ii^ 
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Appendix A: Background field renormalization of 
the Finkelstein term Sf 

In this appendix we present details of the derivation 
of Eq. ((55)) with the help of the background field renor- 
malization. Let us separate the matrix field Q into the 
"fast" (Q) and "slow" (Qo = TV^ATq) modes as 

g ^ Tp-igTo. (Al) 

The effective action for the Qo fields is given by 

exp5eff [Qo] = J V[Q]cxpS[T^^Qn] (A2) 

Since we are interesting in the renormalization of the 
interaction parameters Tab only, we insert the spatial in- 
dependent background field Tq in the action (P^ . The 



result can be written as follows 

SpiT^'Qn] = Sf[Qo] + Sf[Q] + of^'' + of^'' 

.(2)4 , ^(2),2 



where 



(i),i 



+0',>^' + 0't'-^ + Q^, (A3) 
J drJ2 ^abtrl^tabSQtirjabQo, 

<y7i\ab 
an\ab 

JdrY, ^abtrl^^tabQotvA'^^.^.SQ, 

an;ab 

ttT f 

— drJ2 rabtTAZ,a,SQtTA"_^.a,SQ, 

^ J-. 

(A4) 



(2)a 



(2),2 _ 



ttT 



an-.ab 



Or, = inTz / drtiArji 



Here we introduce SQ = Q — A and 

A^^To[7iT„-\ A'^.^a,^ToKtab,T„-\ (A5) 

The effective action 5off[Qo] can be obtained by expan- 
sion of S[Tq^QTq] to the second order in and A^.^^^^ 
Then, we find 

SMQo]-Sf[Qo] = {Of^'') + {Of^'') + ^([OP^'T) 

+i(Or^'^Or)'') + ^([on'> + (0.>, (A6) 



where the average (...) is with respect to action ([3T|) - 
(|32p and we omit terms which do not involve infrared 
divergencies. In general, each term in the right hand 
side of Eq. (jA6p produce contributions which cannot be 
expressed in terms of Qq only. However, all such contri- 
butions cancel in the total expression (jA6[) . Therefore, 
we will not list them below. Expanding SQ in series of 
W according to Eq. (|40p and performing averaging with 
the help of Eq. (|4T|) , we obtain 



ttT f 

SosiQo] = -— drY, Kb^^IntabQt^I^tabQ 

an;ab 

+4:TtTz' j drUriQ (A7) 



where 



iab-^ab + Oiab +('^ab + ab + '^^ ab 



-5T 



(1).2;2 ^ g^r, 



ab 



ab 



(A8) 



and similar for z' . Here the contributions to Fjjj, from 



each term in the right hand side of Eq. (|A6|) are given as 
follows 



cd-.ef 

E^^4"'(-"0, (A9) 
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cd;ef 



(27r) 



(AlO) 



cd;ef 



UJm>0 



(2^ 



(All) 



cci:e/ 



1^ 



l/r^(l).2l2\ ^ ^p(l),2;2 _ 327rr 



and 



cd-.ef 



(27r)2 



[^(cd)^(e/)(^^^) - DDl^''\L,ra)] , (Al3) 



(AM) 



Combing contributions (|M1) - (|A14|) we obtain Eq. ([55]). 

The only non-zero contributions to renormalization of 
z are 



(2^ 



[i^7^p(w„) - (w„)] (A15) 



and 



[Or,) -^5zJ ^ — 2^ Ted 



E ^4^'^^ 



(A16) 



Wm>0 



In total, Eqs. (KT5\ and (|JT6l) give 



647rr 



.5, Ei'^K.). (A17) 



It coincides with Eq. ([51 



Appendix B: Evaluation of DC transconductance a'l^ 



In this appendix we present calculations of the DC 
transconductance in the one-loop approximation. Simi- 
larly to the total conductance, the transconductance can 



be obtained from 

-i2(-«) = (tr[4"t_,Q][/!!„t+,Q]) + ^ 



dr' 



x((tr/:t_0(r)VQ(r)tr/!'„t+Q(r')VQ(r'))) (Bl) 

after the analytic continuation to the real frequencies: 
iiOn Lo + iO^ at (jj — > 0. Here matrices t± = (too i 
t3o)/2. Evaluation of the transconductance according to 
Eq. (|Bip in the one-loop approximation yields 

327rT^ 



ab 



d2p 



(27r)^ 



SnT f d^p 



X I Sp (t+ [tab , t cd] ) t^m EabDp (t^m+n ) £'-0*,°''^ (w,„ ) 

xra;,min{w,„,iLj„ }7^p(cj,„+„)Z?4''^)(w,„) I, (B2) 



where DDp°'''\u!m) = Dp{ujrn)D^"' (uim)- Evaluating the 
traces we find 

a[,{^c.^) = /tI^ E |-™[i?4'"'(^™+n) 

CT^xWn 7 (27r)2 ^j^Q^ L 

-|-min{w,„,i:j„ }p2i^p(c^,„+„)Z?i?^'°'(^m)|. (B3) 

Performing the analytic continuation to the real frequen- 
cies, iujn uj + iO^, one obtains the DC transconduc- 
tance in the one-loop approximation: 

0- 2 = ^Im / 7— / dnl — f^coth — 



Jab) I 



i?420),i?(f^) _p2^2^(20),fl(^)l _^p2Q^^^^ 



2f 



Dl^'^'^'i^)-^ l^Dlm + DZ)(2°)'^(f7)l \, (B4) 



where DDi'''''>'^{n) = D^{n)Di''''^''^{n). Next, Eq. ^ 
can be simplified as 



'12 — 2 



" Re / dr^r^coth — 
2T 



X < z 



n^(i7)l'-(z + r.„)[i?(2")-^^(f}) 



(B5) 
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One can check that due to integration over momentum p ature, a'12 = 0. 
the DC transconductance vanishes at arbitrary temper- 
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